Influence of exciton-exciton correlations on the polarization characteristics of the 
polariton amplification in semiconductor microcavities 
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Based on a microscopic many-particle theory we investigate the influence of excitonic correlations 
on the vectorial polarization state characteristics of the parametric amplification of polaritons in 
semiconductor microcavities. We study a microcavity with perfect in-plane isotropy. A linear 
stability analysis of the cavity polariton dynamics shows that in the co-linear (TE-TE or TM- 
TM) pump-probe polarization state configuration, excitonic correlations diminish the parametric 
scattering process whereas it is enhanced by excitonic correlations in the cross-linear (TE-TM or 
TM-TE) configuration. Without any free parameters, our microscopic theory gives a quantitative 
understanding how many-particle effects can lead to a rotation or change of the outgoing (amplified) 
probe signal's vectorial polarization state relative to the incoming one's. 

PACS numbers: 71.35.-y, 71.36.+C, 42.65.Sf, 42.65.-k 



I. INTRODUCTION 

In the past decade the parametric amplification of po- 
laritons in planar semiconductor microcavities has been 
the subject of intense experimental and theoretical re- 
search, see, e.g., Refs. HlHiili or the reviews given in 
Refs. [7,8. 9. In a typical pump-probe setup in a co-circular 
polarization configuration the amplification of a weak 
probe pulse has mainly been attributed to four-wave mix- 
ing (FWM) processes mediated by the repulsive Coulomb 
interaction of the exciton constituent of the polaritons 
excited on the lower polariton branch (LPB)^^^ For 
a specific pump in-plane momentum (defining the so- 
called "magic angle"), energy and momentum conserva- 
tion is best fulfilled for the FWM processes and thus 
a pronounced angular dependence of this amplification 
is observedi 1 ^ 3 - Since in the strong coupling regime the 
LPB is spectrally well below the two-exciton scattering 
continuum, the influence of excitonic correlations in the 
scattering processes of polaritons on the LPB is strongly 
suppressed [compared to the situation in a single quan- 
tum well (QW) without the strong coupling to a confined 
photon cavity mode 11 ]. However, even for co-circular 
pump-probe excitation these correlations must be con- 
sidered for a complete understanding of the experimental 
results £I2ii2 

Whereas for co-circular pump-probe excitation only 
exciton-exciton scattering in the electron-spin triplet 
channel plays a role, for excitations in other vectorial 
polarization state configurations, excitonic scattering in 
the electron-spin singlet channel is also expected to con- 
tribute to the amplification mechanism. In the latter 
case, a change (in the following loosely referred to as 
'rotation') of the vectorial polarization state of the am- 
plified probe signal compared to the incoming one's can 
be attributed to this coupling of the two spin subsys- 
tems excited with right (+) and left (— ) circularly po- 
larized light, respectively^^^^^^±&2iL2i However, 
different effects can overshadow rotations in the vec- 
torial polarization state that are caused by the spin- 



dependent many-particle interactions that mediate the 
amplification process, e.g., a splitting of the TE and 
TM cavity mode a 22 ' 23 (longitudinal-transverse splitting), 
or an in-plane anisotropy of the embedded QW or the 
cavit y 21 : 24 . Furthermore, for not linearly polarized pump 
excitation an imbalance in the polariton densities in the 
two spin subsystems (+ and — ) will also lead to a rota- 
tion in the vectorial polarization state of the amplified 
signa L 17 : 21 These different mechanisms have previously 
been investigate d 16 : 17 ' 18 : 19 ' 21 based on models describing 
the effective polariton dynamics in the cavity. In these 
models that describe the system dynamics at the po- 
lariton quasi-particle level, the spin-dependent polariton- 
polariton scattering matrix elements are included as in- 
put parameters for the theory. With a reasonable choice 
of the parameter set, good agreement with experimental 
results showing rotations in the probe's vectorial polar- 
ization state has been obtained : 16 ' 18 ' 19 ' 21 

In contrast to these previous studies , 16 ' 17 ' 18 ' 19 ' 21 we 
employ a microscopic theory that calculates, from a 
few material parameters, the scattering matrices driv- 
ing the polariton amplification in the different vectorial 
polarization state channels 2 ^. No additional assump- 
tions for the effective polariton-polariton interaction are 
needed, which is directly included in our theory via the 
frequency dependent and complex exciton-exciton scat- 
tering matrices^ Our theoretical analysis incorporates: 
(i) the well-established microscopic many-particle the- 
ory for the optically-induced QW polarization dynam- 
ics based on the dynamics-controlled truncation (DCT) 
formalis m 26 ' 27 , and (ii) the self-consistent coupling of this 
dynamics to the dynamics of the optical fields in the cav- 
ity mode o 10 ' 28 ' 29 including all vectorial polarization state 
channels. The theory consistently includes all coherent 
third order (x*" 3 "*) nonlinearities and the resulting equa- 
tions of motion are solved in a self-consistent fashion in 
the optical fields which includes a certain class of higher- 
order nonlinearities . 30 ' 31 ' 32 Correlations involving more 
than two excitons and those involving incoherent exci- 
tons are neglected. These effects are not expected to 
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qualitatively alter the presented results for the consid- 
ered coherent exciton densities of ~ 10 10 cm -2 , especially 
for excitation well below the exciton resonance. 

Based on this theory we introduce a linear stability 
analysis (LSA) of the cavity polariton dynamics as a gen- 
eral and powerful tool to study the role of spin-dependent 
polariton-polariton scattering (including time-retarded 
quantum correlations). For steady-state pump excita- 
tion and as long as depletion of the pump from scattering 
into probe and FWM signals can be neglected, the LSA 
gives comprehensive information about growth and/or 
decay (in the following only referred to as 'growth') rates 
for probe and FWM intensities in all vectorial polariza- 
tion states. The growth rates determine the exponential 
growth of components in different vectorial polarization 
states over time, and determine together with the initial 
conditions the ratio of these components after a given 
growth duration, uniquely determining the final vectorial 
polarization state. Although the results are obtained for 
strict steady-state pump excitation, as discussed below 
they can to a large extent be carried over to the analysis 
of pump-probe experiments with finite pulse lengths 

We use this theory to investigate a microcavity sys- 
tem with perfect in-plane isotropy. As an intrinsic ef- 
fect that is not caused by structural imperfection, we 
include a splitting of the TE and TM cavity modes as 
shown in Fig. [TJa). This way we study a system where 
all vectorial polarization state rotations of the amplified 
probe signal can unambiguously be traced back to intrin- 
sic phenomena always present in planar semiconductor 
microcavities, the TE-TM cavity-mode splitting and the 
spin-dependent polariton-polariton scattering mediating 
the amplification. 

For this system we analyze results that show how for 
a linearly polarized pump many-particle correlations and 
the TE-TM cavity-mode splitting lead to different growth 
rates of the linearly polarized components (TE,TM) 
in the probe pulse. This difference can lead to rota- 
tion in the vectorial polarization state of the amplified 
probe compared to the incoming one'si 16 ' 17 ! 18 ' 19 ? 21 Start- 
ing from the equations governing the cavity-polariton dy- 
namics, we take advantage of our theoretical approach to 
isolate and discuss the frequency-dependent scattering 
matrices that give rise to this difference in the vectorial 
polarization state channels. We show that in the studied 
regime, close to the amplification threshold, even where 
the correlation contribution in the spin-singlet channel 
is wea k 15 ! 20 , these correlations can give rise to an al- 
most complete vectorial polarization state rotation into 
the "preferred" cross-linear (TE-TM or TM-TE) pump- 
probe configuration for the pump and outgoing probe 
pulses. For sufficiently long amplification duration, this 
result can become virtually independent of the input 
probe's vectorial polarization state as long as it contains a 
small component polarized perpendicular to the linearly 
polarized pump. 
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FIG. 1: (color online) (a) Shown is the dependence of the 
cavity polariton modes on the magnitude k of the in-plane 
momentum. Depicted are the lower (LPB) and upper (UPB) 
polariton branches for TE (solid) and TM (dashed) cavity 
modes, and the bare cavity and exciton (dotted) dispersions. 
For details on the modeling of the cavity modes see Sec. [TTJ 
(b) Real and imaginary parts of the two-exciton scattering 
matrices T(fi) in the co-circular (H — h) and counter-circular 
(-1 — ) polarization state channels. (a,b) Results are shown for 
typical GaAs parameters^ as used throughout this work. 



II. THE THEORETICAL MODEL 

We use a microscopic many-particle theory to describe 
the coherent QW response to the light field confined in 
the cavity. Based on the dynamics-controlled trunca- 
tion (DCT) approac h 26 ! 27 all coherent optically-induced 
third order nonlinearities, i.e., phase-space-filling (PSF), 
excitonic mean-field (Hartree-Fock) Coulomb interaction 
and two-exciton correlations are included on a micro- 
scopic level. We use a two-band model (including spin- 
degenerate conduction and heavy-hole valence band) to 
describe the optically induced polarization in the GaAs 
QW. 33 Since we are mainly interested in pump excitation 
in the LPB, i.e., energetically below the bare exciton res- 
onance (cf. Fig. [TJ) , we account for the dominant contri- 
butions to the QW response by evaluating the optically 
induced QW polarization in the Is heavy-hole exciton 
basis, 29 ! 31 . 32 ! 34 

We start from the coupled equations of motion for 
the field Ek in the cavity modes with in-plane mo- 
mentum k (treated in quasi- mode approximation 35 ) and 
the optically induced interband polarization amplitude 
Pk in the embedded QW. We formulate our theory in 
the TE-TM basis for the optical fields in the cavity, 
E k = -Bk E e TE + £™e T M (see Fig. [2] for the excitation 
geometry), where the field components in the TE mode 
i?™eTE (also called s-polarized) are characterized by an 
electric field vector with in-plane (in the plane of the 
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plane of 
incidence 




quantum-well 
plane 

FIG. 2: (color online) Schematic of the excitation geometry. 
The plane of incidence is spanned by the wave vector of the 
incoming light field and the z axis. All lines in this plane 
are solid. The quantum-well plane is the x-y plane. All lines 
in this plane are broken. The figure shows the basis vectors 
e TE PlanC alm e TM lane that span the projection of the TE-TM 
basis on the x-y plane. The polar angle d and azimuthal angle 
95k are also shown. For more explanation see text in Sec. [TTJ 



QW) component perpendicular to the in-plane momen- 
tum k, and field components in the TM mode E™eTM 
(also called p-polarized) by an electric field vector with 
in-plane component parallel k. In this basis it is most in- 
tuitive to include different cavity-mode exciton couplings 
for the TE and TM modes in the theory: the in-plane 
component of the fields in the TE mode does not depend 
on the polar angle of incidence t? (cf . Fig. [J) and thus the 
coupling strength to the excitonic dipole in the quantum- 



well plane does not depend on A different result is 
found for fields in the TM mode where the magnitude 
of the in-plane component depends on the polar angle of 
incidence i?. Since the z component of fields in the TM 
mode does not couple to the excitonic dipole for exci- 
tation of heavy-hole excitons in the QW— , the effective 
coupling constant of excitons and fields in the TM mode 
decreases with the polar angle like ~ cosi?. Additionally, 
we include a slightly different polar angular dependence 

TM 

of the bare TE and TM cavity dispersions w fc TE that in 
general depend on the specific materials and design of 
the cavity^ The resulting cavity dispersions are shown 
in Fig.QJa) and the parameters will be given later in this 
section. In order not to complicate the structure of the 
nonlinear terms in the equations of motion for the polar- 
ization amplitudes, we use the usual Cartesian basis (X-Y 
basis) in the QW plane (x — y plane) to decompose the 



polarization into its components as pk = p^e x 



In the X-Y basis, the projection of the TE-TM basis vec- 

■ ,i i in-planc i in-planc 

tors on the x — y plane, 6 TE f and e TIV [ , rotates 
with the in-plane component k of the momentum of the 
incident wave. The azimuthal angle between k and the 
x axis is denoted by iyJk (cf. Fig. \2§ ■ Simple geometric 
considerations lead to the azimuthal angular dependen- 
cies that appear in the terms coupling the equations of 
motion for field and polarization amplitude components 
in the different bases: 
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These equations constitute the generalization of the 
equations given in Ref. 10, now including all vectorial 
polarization states. As discussed above, the polarization 
amplitudes are given in the X-Y basis, while the fields 
are given in the TE-TM basis M- The meaning of the 
symbols in Eq. ([2]) are to be discussed in the remain- 



der of this paragraph, along with the used parameters 
and approximations. Unless otherwise noted, the time 
argument in Eqs. ([1]) and ^ is t. t c is the coupling 
constant of the cavity mode to the external light fields 

TM 

E^f nc , and the dephasing constant j c describes optical 
losses from the cavity to the outside world.— The depen- 
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TM 

dence of the bare cavity modes w fc TE and the dependence 

of the exciton-cavity mode coupling V k TE on the polar 
angle § is modeled on a phenomenological level along 
the guidelines given in Ref. [22. We approximate the bare 
cavity dispersions with w™ = + 100 me V • sin 2 
and = with sintf = ^£2.. This way a TE- 

k cos v ujnbg J 

TM cavity-mode splitting from the mismatch of the 
center of the stopband of the cavity mirrors and the 
Fabry-Perot frequency of the cavity is phenomenologi- 
cally included^ 2 - The cavity-mode exciton couplings are 
V™ = F ™costf and VJ™ = VJ,™, respectively, with 

TM 

V™ = V^ E = 5.2 meV. With fiw TE = eg we assume 
zero cavity-mode exciton splitting for k = 0. The cho- 
sen parameters^ give a reasonable magnitude and polar 
angular dependence of the TE-TM mode split ting in the 
LPB, comparable to the results in, e.g., Ref. |23J. Since 
the presented results do not crucially depend on the de- 
tails of the cavity-mode splitting, no further insight is 
expected from a more elaborate treatment. In Eq. 
Hco^ is the Is heavy-hole exciton in-plane dispersion, ~f x 
a phenomenological dephasing constant of the excitonic 
polarization amplitude, and A is related to the excitonic 
PSF constant A PSF by A = ^gy, with <j>\ s (r) being 
the two-dimensional QW exciton wavefunction. With- 
out loss of generality, in the following the quantities 4>\ a , 

TM 

V TE , and t c are chosen to be real-valued. The parameter 
values are listed in Ref. 33. Although in this paper we 
investigate a spatially isotropic microcavity-system, spa- 
tial anisotropy can easily be included in the theory via 
~ * w k ^0 model an anisotropic dispersion of the QW 

TM TM 

excitons and via w fe TE — > lu^ e to model anisotropy of the 
cavity modes. The two-exciton scattering matrices (T- 
matrices) T in the co-circular (H — h) and counter-circular 
(H — ) polarization state channels include a two-exciton 
dephasing rate 2"f and are given by T ++ = T + and 
T+- = (T+ + T~)/2, with the T-matrices T + and T~ 
in the electron-spin triplet and singlet channel, respec- 
tively, as defined in Eq. (32) of Ref. [25|. The frequency 
dependence of real and imaginary parts of T ++ and T" 1 
is shown in Fig. [TJb). We neglect the momentum depen- 
dence of the T-matrices for scattering processes involv- 
ing two excitons with different in-plane momenta. Cal- 
culations in a different context have shown that this is 
justified in a good approximation for the small optical 
momenta contributing here^ We also neglect possible 
corrections to the excitonic T-matrices from the coupling 
to the photons in the cavity modes. This is supported 
by experimental observations that indicate that even in 
the strong coupling regime the biexciton binding energy 
is not significantly affected by the coupling to the cavity 
modes. 3 - Also, good theory-experiment agreement has 
been achieved in Rcfs. 4fJ using a theory based on 

pure exciton-exciton scattering matrices. 

To simulate a typical pump-probe setup, we start from 
Eqs. (fT]) and ^ and chose a finite in-plane momentum 
k p for the pump propagating along a given axis, here, 



without loss of generality the x axis, i.e., (pk = 0. We 
"detect" the probe in normal incidence with in-plane mo- 
mentum k = where in the past the strongest amplifica- 
tion has been observed) 10 ' 12 ! 16 ' 41 This fixes the in-plane 
momentum of the background- free FWM signal to 2k p . 
We go beyond an evaluation of the theory on a strict 
X^ 3 ) level by self-consistently calculating the resulting ex- 
citonic and biexcitonic polarization amplitude dynamics 
up to arbitrary order in the pump fiel d 30 ' 31 and we lin- 
earize the equations of motion in the weak probe field. 
Only via this self-consistent solution the coupling of the 
probe signal to the background-free FWM signal is in- 
cluded in the theory which provides the basic feedback 
mechanism that leads to the unstable behavior in, e.g., 
Refs. fl.2.3.4.5.7,3- We limit our analysis to coherent ex- 
citon densities that are low enough [<2x 10 10 cm -2 , cf. 
Figs. [3] and [4] so that neglect of higher than two-exciton 
Coulomb correlations can be justified. 



III. LINEAR STABILITY ANALYSIS 

In this section we introduce the linear stability analysis 
(LSA) used in the remainder of this paper. To analyze 
the stability of the pump-probe dynamics the LSA is done 
without an incoming probe field and for a linearly po- 
larized monochromatic continuous wave (cw) pump field 

TM _ TM 

^k TB (*) = ^k TE {^p)G~ l ^ pt inducing the pump polariza- 

X X , .X 

tion amplitude (t) — p~7 (ui p )e~ iu ' pt , with p7 = and 

X TM 

E-^ E — (u>p is the pump frequency). The pump po- 
larization amplitude is a solution of the cubic nonlinear 
pump equation following from Eqs. (fTJ) and ([2]) for unidi- 
rectional light propagation and is determined as outlined 
in Appendix |Al The resulting coherent exciton density 

X 

\p£ I for excitation with a linearly polarized pump of 
fixed intensity is shown in Figs. [3] and Has a function 
of the magnitude of the pump in-plane momentum |k p | 
and pump detuning Ae from the bare exciton resonance 
for excitation in the TE or TM mode, respectively. No 
bistable behavior of the pump-induced exciton density 
in Figs. [3] and [4] is found (this follows from the solu- 
tion of the nonlinear pump equation as outlined in Ap- 
pendix[A| • However, for other values of cavity or QW pa- 
rameters and a different pump intensity, bistability may 
occur, which would complicate our discussion 41,42 . The 
pump densities shown have their maxima close to the lin- 
ear polariton dispersions (included as the dashed lines) 
and decrease along the LPB with increasing in-plane mo- 
mentum because of the decreased coupling of the (mostly 
exciton-like) large-momentum polariton states to the in- 
coming field. Furthermore, a reduced exciton density 
is found for excitation of the UPB caused by strong 
excitation-induced dephasing (EID) for pump excitation 
in the two-exciton scattering continuum (the spectral re- 
gion where the two-exciton scattering matrices shown in 
Fig. [U exhibit a large imaginary part). For the stability 
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analysis we evaluate the memory integrals in Eq. ((2|) con- 
tributing to the probe and FWM directions in a Markov 
approximation for the two-exciton scattering continua in 
the T-matrices T ++ and 7~ H . Our Markov approxi- 
mation is effected by taking po(t') « po(*)e I " p< -'~ t ' and 
P2k p (t') ~ P2k p (t)e iuJp( - t ~ t ' ) where the probe and FWM 
polarization amplitudes p and P2k p appear under the 
time-retarded integrals in Eq. @ together with the con- 
tinuum contributions in T ++ and . In contrast, we 
include the bound biexciton state exactly via the time- 
dependent amplitudes bk (t) and 63k p (t)- For this we 
separate the bound biexciton contributions from the 
correlation kernels T+~ in Eq. © as T+~ = T+~ t +T+T. 
The bound biexciton contributions to Eq. |2]) can be ex- 
actly included via the equations of motion [cf. Eqs. (11) 
and (12) of Ref. l25[ for the biexciton amplitudes bk p (t) 
and &3k (t) which are labeled according to the total in- 
plane momentum of their source terms: ~ poPk p and 
~ P2k p Pk p , respectively. This way we include quantum 
memory effects related to the excitation of bound biexci- 
tons, which were previously shown to play an important 
role in the study of FWM instabilities in single semicon- 
ductor QWsii. Since, for the chosen excitation geometry 
with finite pump in-plane momentum k p , the probe and 
FWM signals do not oscillate at the pump frequency uo p , 
the Markov approximation for the two-exciton scattering 
continuum may not be as justified as it is for the single 
QW system investigated in Ref. [ll|. However, close to or 
in the unstable regime, those wave mixing processes that 
describe the pairwise scattering of pump polaritons into 
the probe and FWM directions play the dominant role in 
the probe and FWM dynamics. For monochromatic cw 
pump excitation these terms are of purely Markovian na- 
ture and hence the T-matrices in these terms contribute 
exactly at frequency Q — 2uj p ; no approximation is re- 
quired. And indeed, for the results discussed here, not 
even from the excitation of the bound biexciton state 
have we found a sizable contribution from non-Markovian 
(quantum memory) effects. For monochromatic cw pump 

TM „ TM 

excitation and with the ansatz E^ E (t) = -E TE (i)e -Mtf » ,t , 

TM ~ TM XX X 

^ 2 T k E p (*) = £ 2 T k E p (*)e-^*, P J (*) - P Y (t)e— V, p- K (t) = 
fipWe-^' and b kp (t) = b kp (t)e^ pt , b 3kp (t) = 

b3k p (t)e~ t2uJpt the coupled probe and FWM dynamics can 
be written in the form 

Hp(t) = Mp(t) . (3) 



The vector 

p(t) = [E™(t),E™;(t),E™(t),E™*(t), 

PoW^tkpW^^W^^W^kpW^kpW] 7 , 

groups field and polarization amplitude variables to- 
gether. M is a time-independent matrix where all system 
parameters and the steady-state pump polarization am- 
plitude (Figs. [3] and [4]) and the corresponding pump field 
in the cavity modes parametrically enter the analysis. 
For excitation with a linearly polarized pump either 
exciting the TE or the TM mode, the matrix M is block- 
diagonal for the components of probe and FWM parallel 
(co-linear configuration) and perpendicular (cross-linear 
configuration) to the pump's vectorial polarization state. 
Then for each pump polarization state (TE or TM, 
respectively), the 10 x 10 matrix M can be decomposed 
into the 6x6 block M^ s ' <Tp with a s = a p and the 4x4 

block Mj_ 3,<Jp with a s 7^ <j p which describe the dynam- 
ics of the coupled variables in the vectors p,f s (t) = 

IK' (t),E££(t),pg- (t),p£* p (t),b k M k K (t)] T 

with a s = a p for the co-linear configurations and 
Pl'{t) = {K°(t),E^(t),p^(t),p£*(t)} T with a s ± a p 
for the cross-linear configurations, respectively. The 
indices a s and a p relate to the cavity modes TE and 
TM, or to the corresponding excitonic polarization 
amplitude components X and Y that are excited by the 
fields in these modes. Note, that for the above-described 
excitation situation (pump pulse propagating along the 
x axis and linearly polarized excitation in TE or TM 
mode) the x component of the polarization amplitude 
is exclusively excited by fields in the TM mode and the 
y component by fields in the TE mode. The matrices 
M|j Ts,crp and M° s ' <Jp can be derived from Eqs. fl} and 
@ and take the following form: 
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The time-independent coefficients are defined as: 
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T ++ (il) and T + (f2) are the Fourier transformed cor- 
relation kernels as shown in Fig. [TJb) and defined by 
Eqs. (27) and (32) of Ref. U The coupling strength 
C xx w 0.54 Efa^ of the excitonic polarization amplitudes 
to the bound biexciton amplitude is given by the biexci- 
ton ground state wave function £(q) in the electron spin 
singlet configuration and the corresponding two-exciton 
Coulomb interaction matrix element W~ x (q,0), both as 
defined in Eqs. (14) and (24) of Ref. [U 

For steady-state monochromatic pump excitation, the 
linear stability analysis formulated in this section gives 
us comprehensive information about growth (real-part 
of the eigenvalues of M) and frequency (imaginary-part 
of the eigenvalues of M) of the polariton modes in the 
probe and FWM directions, after the initial external driv- 
ing pulse (seed) in the probe direction is gone. Infor- 
mation how the different eigenmodes contribute to the 
polarization amplitudes and fields with different in-plane 
momenta (0 or 2k p ) can be obtained from the eigenvec- 
tors of the matrices M, If at least one of the eigenvalues 
Aj of the matrices M fulfills Re{Ai} > 0, the system is 

X X 

unstable. An arbitrarily small seed of p$ or p^. (X or 
Y, depending on which subspace shows an unstable dy- 
namics) or in the corresponding cavity modes would grow 
exponentially until the matrix M ceases to describe the 
system correctly. 

In addition to the strict steady-state analysis (regard- 
ing the pump excitation), the general information ob- 
tained from the stability analysis can - to a large extent 
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FIG. 3: (color online) Top: Coherent steady-state exciton 
density IPfcJ 2 for a fixed pump intensity of a linearly TE po- 
larized pump as a function of the magnitude k p of the pump 
in-plane momentum and the pump detuning Ae = hui p — eg 
from the bare exciton resonance eg. The linear polariton dis- 
persions are included as the dashed lines. Middle and lower 

TE TE 

figures show the real part of the eigenvalue of My ' and 

TE TIvI 

M ± ' with the largest real part in meV (maximum growth 
rate if larger than zero). The linear polariton dispersions are 
included as the dashed lines and the insets show the same 
data around the "magic angle" for pump excitation on the 
LPB. 



- be carried over to the discussion of pump-probe experi- 
ments with finite pulse lengths. From the linear stability 
analysis reasonable predictions for growth rates of probe 
and FWM signals and thus for polarization rotations can 
be made as long as no external probe pulse significantly 
drives the probe polariton dynamics during the period of 



7 



coherent steady-state exciton density [10 1 °cm 2 ] 



E 





■ ' J> 




5- 


TM 









0.3 






0.5 






9? 






1.1 
1 3 


0- 


TM ■ 


— 1.6 

1.8 






^™ 2.1 


-5- 
t 






11 2 3 4 5 





pump in-plane momentum k [10 6 m 1 ] 




2 4 6 

pump in-plane momentum k [10 6 m 1 ] 



FIG. 4: (color online) Same as Fig. [3] but for linearly TM 
polarized pump. 



amplification. Furthermore, for interpretation of pulsed 
experiments based on the linear stability analysis, the 
pump pulse must be spectrally sufficiently narrow and 
pump and probe must have significant temporal overlap. 
After a sufficiently long period of time that particular 
eigenmode of M corresponding to the eigenvalue with 
the largest real part will dominate the overall outgoing 
signal in probe and FWM directions. In a pump-probe 
experiment an incoming probe pulse in this particular 
mode will be most efficiently amplified, or for steady- 
state pump excitation without an incoming probe, fluc- 
tuations in this mode (serving as a seed) will grow most 
efficiently over time and dominate the signal in probe and 
FWM direction after a sufficiently long growth period. 
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pump in-plane momentum k [10 6 m 1 ] 



FIG. 5: (color online) Same as the two lower panels in Fig. [3] 
but without exciton-exciton scattering in the +— channel 
(f +" ee 0). 



IV. RESULTS AND DISCUSSION 

Without excitonic correlations (and neglecting the TE- 
TM-splitting of the cavity modes) the equations of mo- 
tion for the two circular polarization state channels, + 
and — , are decoupled. In this case, for excitation with a 
linearly polarized pump, where equal densities of polari- 
tons are excited in these two different polarization state 
channels, the incoming (seed in our linear stability anal- 
ysis) and outgoing probe are always in the same vectorial 
polarization state, in the stable (all Re{A} < 0) as well as 
in the unstable (at least one Re{A} > 0) regime. Neglect- 
ing the TE-TM cavity-mode splitting, for the pump in 
a linear polarization state, only the spin-dependent exci- 
tonic correlations can give rise to a rotation of the outgo- 
ing probe signal's vectorial polarization state relative to 
the incoming one's. The actual fraction of polaritons that 
is scattered into the probe and FWM directions in a vec- 
torial polarization state perpendicular or parallel to the 
pump's vectorial polarization state, respectively, strongly 
depends on the excitonic correlations in the H — polar- 
ization state channel (included in T H ).= As discussed in 
the previous section, for a linearly polarized (either TE or 
TM) pump, all the eigenmodes [eigenvectors of the ma- 
trix M in Eq. ([3])] for the probe and FWM dynamics are 
either polarized parallel (co-linear configuration) or per- 
pendicular (cross-linear configuration) to the pump, even 
when both excitonic correlations and cavity-mode split- 
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ting are included. However, either excitonic correlations 
alone or the cavity-mode splitting alone (when the po- 
lariton scattering is mediated by Hartree-Fock Coulomb 
interaction) is sufficient to give different probe and FWM 
dynamics in the two different polarization state configu- 
rations. 

Since the above-listed effects lead to a difference in the 
growth rates of the modes polarized parallel or perpen- 
dicular to the pump, for an arbitrarily polarized probe 
the two different vectorial polarization state components 
(parallel or perpendicular) will grow differently over time. 
The stronger the amplification of the probe and the 
longer the amplification duration, the more the fractions 
of the probe in those modes that exhibit the fastest ex- 
ponential growth, will dominate the outgoing probe and 
FWM signals. Thus, for strong amplification, the growth 
rates of the fastest growing modes in the two polarization 
states (parallel and perpendicular) ultimately determine 
the rotation in the probe's vectorial polarization state. 
Studying these growth rates also answers the question 
about the preferred mode for the growth of probe fluc- 
tuations over time, when no incoming probe is present. 
Figures [3J and 0] show these growth rates - the real part 
of the eigenvalue of M with the largest real part in each 
case - for a fixed pump intensity for the different po- 
larization state configurations (co-linear or cross-linear) 
for pump excitation of TE (Fig. \5$ and TM (Fig. EJ 
mode, respectively. The results show that pumping ei- 
ther TE or TM mode does not significantly influence the 
overall result regarding the maximum growth rates in 
the co-linear or cross-linear configuration; merely a small 
change in the optimum pump momentum and frequency 
is observed. However, a significant difference between co- 
linear and cross-linear configuration is found for pump- 
ing close to the inflection point of the LPB (the so-called 
"magic angle" ) where phase- matching is best fulfilled so 
that triply-resonant (resonant for the pump excitation 
and at an angle so that the dispersions in probe and 
FWM directions allow for phase-matched scattering of 
pairs of pump-excited polaritons into these two direc- 
tions) amplification of the polaritons can occur. For the 
chosen intensity, close to the instability threshold, and 
for pump excitation under the "magic angle" we are in 
a regime where instability (Re{A} > 0) and correspond- 
ing exponential signal growth is only found in the cross- 
linear configurations (TE-TM and TM-TE) while all the 
modes in the co-linear configurations (TE-TE and TM- 
TM) are exponentially decaying. In this regime, for an 
arbitrarily polarized probe, only that component polar- 
ized perpendicular to the pump is exponentially grow- 
ing over time and thus only this component experiences 
significant amplification. To isolate the mechanism that 
leads to this striking difference in the polarization state 
configurations, Fig. [5] shows results for both configura- 
tions in Fig. [3J but without taking into account the cor- 
relations in the H — channel (T H = 0). Without these 
correlations the results for the two configurations almost 
look alike; we find only a small difference in the growth 



rates caused by the TE-TM cavity-mode splitting. Note 
that without the correlations in the H — channel in both 
configurations the instability threshold is not reached for 
the same pump intensity as used in Figs. [3J and [3J 

In the co-circular (-1 — (-) excitation configuration it was 
found earlier— that excitonic correlations in the ++ 
channel may considerably reduce the maximum growth 
rates in the polariton amplification. In the H — h channel 
the driving mechanism for the instabilities, the phase- 
conjugate feedback, is weakened by the two-exciton cor- 
relations. Additionally, but for large negative detuning 
less important, correlations in the H — h channel give rise to 
pump- induced EID, also reducing the exponential growth 
rate over time. Figures [3J and E] show that the correla- 
tions in the H — channel enhance the growth rate in the 
cross-linear configuration and diminish it in the co- linear 
configuration, compared to the results shown in Fig. [5] 
where these correlations are absent. 

In the following we will interpret these results in 
terms of the exciton-exciton scattering matrices shown 
in Fig. [TJb) for the different polarization state channels 
of polariton-polariton scattering. For this discussion we 
ignore the small PSF nonlinearities that contribute to 
the probe and FWM dynamics, and concentrate on the 
nonlinearities in the excitonic polarization amplitudes in 
Eq. @ that contribute to the probe and FWM dynamics 
and therefore determine the amplification process. Being 
sufficient for a qualitative understanding we discuss all 
contributions in Markov approximation. Three different 
terms have to be analyzed that enter the matrices M in 
Eqs. gl) and ©: 

(i) Excitation-induced dephasing for the excitonic com- 
ponent of the polaritons (entering M via CLuV)' 

Im{f ++ (2a, p ) + ^.||t+-(2c p )}|^| 2 . (6) 

(ii) Nonlinear shifts to the effective exciton resonances in 
probe and FWM direction (entering M via a^): 

Rc{r++(2^) + 4||r+-(2^)}|^| 2 . (7) 

(iii) The phase-conjugate oscillation feedback for the ex- 
citonic constituents of the polaritons that drives the in- 
stability (entering M via b" p ): 

i((-l)<^f++(2u; p ) +f + -(2u p ))pf . (8) 

The index i e {_L, ||} labels the co- and cross-linear polar- 
ization state configurations, respectively. The maximum 
growth rates in Figs.[3jto[5]are obtained when the pump is 
tuned about 3 meV below the exciton resonance and close 
to the "magic angle". For this pump detuning and in 
Markov approximation the two-exciton scattering matri- 
ces shown in Fig. [TJb) contribute at Kl — 2e x « —6 meV. 
For this detuning the Im{T ++ } is much smaller than 
lm{T' 1 } which according to Eq. ^ leads to a much 
larger EID in the co- linear (||) configuration. According 
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to Eq. ([7|) a partial cancelation of the nonlinear energy 
shifts from contributions in the H — h and H — channels is 
found in the co-linear (||) configuration only. This, how- 
ever, only slightly modifies the effective resonance fre- 
quencies (polariton dispersions) and thus slightly changes 
the optimum pump momentum and frequency. Most im- 
portant for the stimulated amplification process is the 
difference in the polarization state configurations that 
can be seen in Eq. (|8j). Whereas for the co- linear (|j) 
configuration the sum of the scattering matrices in the 
H — h and H — polarization state channels determines the 
strength of the phase-conjugate feedback driving the in- 
stability, for the cross-linear (_L) configuration the dif- 
ference of the two is relevant. Compared to the results 
without correlations in the H — channel, the difference 
in sign in the real parts of these two contributions [cf. 
Fig. QJb)] leads to strong cancelation in Eq. ([8]) for the 
co-linear configuration and to strong enhancement in the 
cross-linear configuration. As previously pointed out in 
Ref. [ID this leads to an imbalance in the pairwise scat- 
tering of polaritons into the probe and FWM directions 
with polarization parallel or perpendicular to the pump. 
From our results we conclude that this does not necessar- 
ily lead to an overall rotation of the vectorial polarization 
state by 90° in contrast to the conclusions in Ref. [2ll . 
Based on our explicit treatment of the bound biexciton, 
we also find that in the studied system the scattering 
matrix element that drives the amplification pro- 

cess by scattering of polaritons with opposite spins is not 
small compared to 7~ ++ . We note that this is in contrast 
to the case reported in Ref. [H. Our results indicate that 
in the cw regime for a linearly polarized pump and close 
to threshold, spontaneous fluctuations preferably grow 
in the cross-linear polarization configuration. This is in 
agreement with recent observations for a slightly different 
system and excitation geometry.— 

The above discussion has a very general character. 
The actual cancelation of the different contributions in 
Eqs. I©-© depends quantitatively on parameters such 
as the coupling strength of the cavity-modes to the exci- 
tonic polarization amplitudes. The general trend follows 
from the frequency dependence of the two-exciton scat- 
tering matrices as shown in Fig. [ljb): The stronger the 
cavity-mode exciton coupling (shifting the "magic angle" 
to larger negative detuning) the less pronounced the role 
of T H will be. However, especially close to threshold 
even a small difference in growth rates can be crucial, and 
in this regime even a small contribution can play a 
major role for the analysis of vectorial polarization state 
rotations. Although for the system studied here it was 
found to be almost insignificant, the overall role of the 
TE-TM cavity-mode splitting can depend on parameters 
and excitation conditions, too. Furthermore the impor- 
tance of a TE-TM cavity-mode splitting can be different 
in other systems such as the quasi-one-dimensional mi- 
crocavity system studied in Ref. [20L 

Although not relevant for the amplification of polari- 
tons on the LPB, we finally note that Figs. [3] to [5] show 



strong pump- induced E1D for excitation on the UPB. Al- 
though the two-exciton scattering matrices evaluated in 
the Is approximation are not quantitatively accurate in 
this energy region, even in a theory including only coher- 
ent excitations this additional EID would likely inhibit 
the observation of any instability in the UPB (previously 
discussed, e.g., in Ref. [43h in analogy to the situation for 
positive pump detuning in a single QW 11 . 



V. CONCLUSIONS 

Based on a microscopic many-particle theory we have 
investigated the influence of excitonic correlations on the 
vectorial polarization state characteristics of the para- 
metric amplification of polaritons in semiconductor mi- 
crocavities. By means of a linear stability analysis it 
has been analyzed how a linearly polarized pump can in- 
duce a polarization state anisotropy in an otherwise per- 
fectly isotropic microcavity system. For the discussion 
of this effect we take advantage of our theoretical ap- 
proach which - in contrast to previous model o 16 ! 17 ' 18 ? 19 
- is based on microscopically calculated exciton-exciton 
scattering matrix elements. Accounting for all coherent 
correlations between two excitons, these matrix elements 
determine the nonlinear cavity-polariton dynamics in the 
probe and FWM directions in the amplification regime. 

A previous study-i^ found that excitonic correlations 
weaken the polariton amplification for co-circular pump- 
probe excitation. We confirm these findings and addi- 
tionally investigate the effects of correlations on the po- 
lariton amplification in linearly polarized pump-probe 
configurations. We find that scattering contributions 
of the excitonic components of polaritons with opposite 
spins can strongly diminish the driving force for the am- 
plification, the phase-conjugate coupling, in the co-linear 
(TE-TE or TM-TM) pump-probe polarization configura- 
tion and strongly enhance it in the cross-linear (TE-TM 
or TM-TE) configuration. In the spectral region where 
instability occurs the scattering of polaritons with oppo- 
site spins is dominated by the virtual formation of bound 
biexcitons. 

In or close to the unstable regime this polarization 
state anisotropy has the potential to alter the polariza- 
tion state of an amplified probe pulse compared to the in- 
cident probe. If the incoming probe has both components 
polarized parallel and perpendicular to the pump's po- 
larization state, in general the maximum growth rate for 
these two components over time is not the same. Then af- 
ter a certain growth period the amplification of these two 
polarization state components will be different, and hence 
the vectorial polarization state of the outgoing probe is 
rotated compared to the incoming one's. Since the probe 
component polarized perpendicular to the pump's polar- 
ization vector grows faster over time than the parallel 
component, the probe is always rotated toward the "pre- 
ferred" cross-linear configuration. However, the overall 
rotation in the vectorial polarization state depends on 
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both the duration of amplification and the difference in 
the growth rates in the two polarization state channels. 

Finally, we note that in a situation where steady-state 
pump excitation brings only the cross-linear configura- 
tion above the unstable amplification threshold, without 
an incoming probe, only spontaneous fluctuations in the 
cross-linear polarization state channel will be amplified 
and thus observed as a finite signal in probe and FWM 
direction. 
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APPENDIX A: NONLINEAR PUMP EQUATION 

The linear stability analysis of probe and FWM dy- 
namics in Sec. Mil is done for monochromatic cw pump 
excitation. The stationary pump field inside the cav- 

TM 

ity, E^ E , and pump-induced polarization amplitude in 

the QW, p£ , which enter the matrix M in Eq. are 
needed to analyze the probe and FWM dynamics. In this 
work we have considered pump excitation with a linearly 
polarized pump in the TM or TE cavity mode with in- 
plane momentum k p along a certain axis, here, without 
loss of generality the x axis, i.e., ipk p — 0. Seeking a 
steady-state solution for the pump-induced polarization 



Eq. 



TM 
TE 



~ TM 

TE ^ — ILUjjt 



E. 1 ". e 

k D ,mc 



amplitude p£ we use the ansatz E£ 

TM „ TM 

E^ E = E^ E e~ luJpt for the incoming pump field and 
the field in the excited cavity mode, respectively, and 

XX. 

Pk ~ Pk e~ tuJpt for the polarization amplitude, with 

^ TM ^ TM . X 

^ViBc = E ^ = K p = °- Thcn from E q- © it: fol " 

lows that the electric field in the cavity mode for in-plane 
momentum k p is given by: 



~ TM 

= 



TM X _ TM 



tlllJr, 



TM 



(Al) 



For unidirectional light propagation, i.e., by removing all 
sums in Eq. taking all field and polarization ampli- 
tude variables at momentum k p , and replacing the pump 
field by Eq. (jAlj) , a cubic equation of the form 



= a + ax \p£ 



I 2 | 1 4 I x 

I + a 2\Pk v \ +\pL 



(A2) 



can be derived. This equation determines the monochro- 
matic solutions for each pump frequency uj p and incom- 

TM 2 

ing pump intensity (~ \E™.\ ). 



arc 



a = 



ITPV 



J. 1V1 — «w _L IV* 

fc2j-2T/TE PTE 
1 11 l C V k„ l^k„,i] 



\T\- 



TM 



) 2 +7c 2 



TM 2 ^ _ TM 

2hHlV k : E A\E™.f 



and 



«2 



(Re{e}Re{f } + Im{e}Im{f } 



TM 2 _ „ TM 

hHlV™ A2\E- E in f 



with the definitions 



TM 2 

VJ E 



TM 

Hujp — fuVu E + ^7c 



T 4 



TM 

hujp — Hu)^ E + i"f c 



The coefficients in 



Being a cubic equation in \p£ | 2 , depending on the coef- 
ficients ao,ai,a,2, Eq. (|A2)> can have either one or three 
real- valued solutions. These solutions are given, e.g., in 
Ref. 0. 

Equation (|A2|) only determines the magnitude of the 

X 

pump-induced polarization amplitude \p£ | and contains 
no information about the different phases of the polariza- 
tion amplitude and the incoming field. We define a phase 

„ TM X _ TM X 

according to Eg^-p^ = \E^J\p^\e^. This phase 
<j) is required to determine the field in the cavity mode 
from Eq. (|A1[) that is coupled to the polarization ampli- 
tude. The field enters the linear stability analysis via the 
PSF terms in M. If we choose the solution of Eq. (|A2[) 

X X 

to be real, p£ = \p£ |, then the incoming field inducing 

~ TM _ TM 

this polarization amplitude is E^ E jnc = |^k rB inc |e l( *'. The 
phase of the incoming field can be obtained from Eq. ([2|) 
and is given by 

X X TM 

# _ (fig, I + T K\ 3 ) - + !2j 

^ TM X . . ~ TM 5 

ifceVg* (i-^Kj 2 )|£; k T p E inc | 

which then also determines the phase of the field in the 
cavity mode given by Eq. {MJ. 
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